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Abstract 

Stacks of a few intrinsic tunnel junctions were micro- fabricated on the sur- 
face of Bi2Sr2CaCu208+5 single crystals. The number of junctions in a stack 
was tailored by progressively increasing the height of the stack by ion-beam 
etching, while its tunneling characteristics were measured in-situ in a vac- 
uum chamber for temperatures down to ~13 K. Using this in-situ etch- 
ing/measurements technique in a single piece of crystal, we systematically 
excluded any spurious effects arising from variations in the junction parame- 
ters and made clear analysis on the following properties of the surface and in- 
ner conducting planes. First, the tunneling resistance and the current-voltage 
curves are scaled by the surface junction resistance. Second, we confirm that 
the reduction in both the gap and the superconducting transition temperature 
of the surface conducting plane in contact with a normal metal is not caused 
by the variation in the doping level, but is caused by the proximity contact. 
Finally, the main feature of a junction is not affected by the presence of other 
junctions in a stack in a low bias region. 



1 



I. INTRODUCTION 



Recently both superconducting and normal-state tunneling characteristics of junc- 
tions intrinsically formed in crystals of extremely anisotropic high-Tc materials such as 
Bi2Sr2CaCu208+5 (Bi-2212) and Tl2Ba2Ca2Cu30io+5 (Tl-2223) have attracted much re- 
search attention. Previous experimental investigations by many groups have revealed 
that the superconducting order parameter in the crystals is periodically modulated along 
the c-axis and the interlayer coupling at low-enough temperatures has a Josephson nature. 
Existence of such naturally formed intrinsic Josephson junctions has been directly confirmed 
by c-axis tunneling measurements, using small-sized stacks or mesa structures on the sur- 
face of single crystals, including current-voltage {IV) characteristics. 

Hi 

temperature or 

magnetic-field dependence of the critical current ,EHli dynamics of Josephson vortices,EHHi 
microwave responses and the Josephson microwave emission.§ 

Electron tunneling spectroscopy has been conceived as one of the most powerful means 
to determine superconducting gap and its anisotropy.i2hil Compared to other spectroscopic 
techniques,li3^ electron tunneling spectroscopy has a very high energy resolution of the 
electronic structure near the Fermi level Ep. Tunneling spectroscopy and its interpretation 
on high-Tc materials using conventional techniques such as scanning tunneling spectroscopy, 
however, have often been controversial, due to problems of poor surface characterizations 
mainly arising from an extremely short coherence length and uncertainty in the transverse 
momentum direction of the particles upon tunneling via a small t ip. ilil These problems 
have been conveniently circumvented by a tunneling study using intrinsic junctions (Us) 
in high-Tc materials. Especially, tunneling measurements on a small-sized stack, including 
only a few intrinsic junctions, provide valuable information on the nature of the inherent 
superconducting gap and the pairing mechanism of unusual symmetry in the materials. 

It is unsettled yet, however, whether the individual junctions in a stack behave inde- 
pendently or are correlated with each other. Junctions in a stack could be inductively 
coupled by supercurrents flowing on thin conducting planeJi^E or coupled by charging of 
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thin conducting planes by tenneling electrons.Ei Thus the plasma oscillations excited in a 
stack of intrinsic junctions, for example, show collective modes-IHl However, some previous 
studies indicate the opposite may be true. For dc bias current applied along the c-axis, 
coupling between adjacent junctions is small enough so that each junction can be assumed 
to be independent .Si Correlated junction behavior would affect the usefulness of a stack for 
device applications. For instance, correlated junction behavior is useful for high-frequency 
oscillator applications but detrimental to voltage standards applications.il Thus, the issue 
of correlation of the tunneling characteristics of intrinsic junctions still remains of prime 
concern from both academic and device-application points of view. 

Tunneling spectroscopy on the surface of high-T^ materials is widely used to probe the 
nature of order-parameter symmetry, mostly using hybrid junctions consisting of conven- 
tional and high- Tc superconducting electrodes.SB In this case, it is utterly important to 
make sure that the superconducting properties of the surface layer, such as the transi- 
tion temperature, the gap value, and its order-parameter symmetry, do not deviate from 
those of inherent bulk properties.^ Recently, surface-sensitive spectroscopy such as anlge- 
resolved photoemission spectroscopy,ii scanning tunneling spectroscopy,0ii break-junction 
tunneling ,@ and intrinsic-junction tunneling is increasingly used to investigate the 
in-plane normal-state conducting state of high-Tc materials i.e., the existence and nature 
of the pseudogap, especially as varying the doping level. Since the size of the pseudogap is 
known to be sensitive to the doping level in the underdoped regime,&@ any possible sig- 
nificant oxygen deficiency in the surface layer would lead to a very serious misinterpretation 
of the intrinsic bulk properties of the conducting planes. Additionally, our previous study0 
indicates that depositing noble-metal film such as Au on the surface of high-Tc materials, 
often used for surface protection for electron conduction measurements, causes a suppression 
of the surface layer's superconductivity by its proximity contact to normal metal. In this 
respect, it is of prime importance to probe any variation in the doping level of the surface 
layer from that of the inner layers with intrinsic bulk properties. 

Most studies on the intrinsic tunneling effect in high-Tc superconductors, both in super- 



conducting and normal states, have been done to date by fabricating a small elevated stack 
structure (with the cross-sectional area usually of a few hundreds /im^) on the surface of 
single crystals with a few IJs in it. The number of the junctions in a stack is varied by 
changing the height of the stack during the fabrication process using ion-beam etching or 
chemical wet etching. In this way, however, accurate tailoring of the number of junctions 
is a difficult task. To get deeper understanding of the interlayer coupling characteristics 
of an individual intrinsic junction as well as the conducting properties in a plane, both in 
superconducting and normal states, careful in-situ control of the number of IJs in a stack is 
required while monitoring their tunneling characteristics.! 

In this report, we first examine the possibility of any significant variation of the doping 
level in the surface layer of high- Tc single crystals in contact with normal metals, along with 
a change in the superconducting properties. We then discuss the physical implication of the 
scahng of tunneling resistance and the IV curves with respect to surface junction resistance. 
We also examine whether the main feature of a junction is infiuenced by the presence of other 
junctions in a stack. To that end, we focus on the evolution of the tunneling characteristics 
of IJs while progressively increasing the number of junctions in a stack fabricated on the 
surface of Bi-2212 single crystals. Ion-beam etching was used in combination with in-situ 
cryogenic measurements, similar to the technique reported in Ref. |. The number of IJs in 
a stack was increased in sequence by carefully controlling the low-energy ion-beam etching 
time with preformed metallic electrodes on the surface of the stack. Between each etching 
stage the c-axis tunneling IV characteristics and the temperature dependence of the c-axis 
tunneling resistance Rc{T) were taken while varying the substrate temperature from room 
temperature down to ~ 13 K. 

This study directly confirms our previous finding that a weak junction (WJ) forms at 
the surface of a Bi-2212 stack0 which is in contact with a normal-metal electrode. The 
in-situ etching/measurements technique confirms that the suppression of the gap and the 
transition temperature in the surface layer of a stack does not result from the change in 
the doping level or from the proximity- induced weak superconductivity in the Bi-0 layeri'i 
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which may be surfaced. Rather, the suppression results presumably from the proximity 
contact of the surface layer with d^2_y2 symmetry to the Au normal metal. This study also 
convincingly demonstrates that the individual junction is little influenced by the presence of 
other junctions in a stack in a low bias region. As to the cause of development of negative 
dynamic resistance in the high-bias region of the IV curves with increasing number of 
junctions in a stack, we find that the Joule heating effect, not the nonequilibrium, can be a 
good candidate for that cause. 

II. EXPERIMENTS 

Bi-2212 crystals were grown from a melt using alumina crucibles, details of which are 
described elsewhere.0 Pieces of single-crystal platelets with typical size ~ 0.8 x 0.4 x 0.03 
mm^ were selected from a mass of the cooled melt. A platelet was glued onto a sapphire sub- 
strate by negative photoresist, then hard-baked. An optically smooth surface was prepared 
by cleaving a platelet of single crystal using Scotch tape. Right after cleaving, 1000-A-thick 
Au film was deposited on the surface of the Bi-2212 crystals to protect the surface from any 
contamination during stack-fabrication processes. 

We then patterned a large 'base mesa' of size ~ 450 x 15 x 1 fim^ on the crystal surface, 
using standard photolithography and Ar-ion-beam etching, with beam voltage and current 
density Vbeam = 300 V and heam = 0.8 mA/cm^, respectively. Any residue of burned- 
out photoresist on the crystal surface was stripped off by oxygen-plasma etching. Then a 
patterned layer of photoresist was placed around the base mesa and hard-baked to insulate 
the region, excluding the top surface of the base mesa. A thick Au film (~ 8000 A) was 
deposited,^ and connection pads from the base mesa were patterned and ion-beam etched. 
Thereafter, the electric pads of Au (totally ~ 9000 A) crossing the base mesa acted as masks 
for the small stacks in the further in-situ etching process. In this stage, a thin (less than 300 
A in nominal thickness) Au film was intentionally left on the top surface of the base mesa 
between the pads to prevent the inadvertent premature formation of small stack structures. 
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The remaining Au film would be removed during the next in-situ measurement process with 
controlled ion-beam etching. Figure 1 illustrates the optical photograph of a specimen after 
completing the whole processes prior to the in-situ etching. 

For in-situ monitoring of the tunneling characteristics of a small stack, the specimen 
was loaded into an ultra-high- vacuum chamber equipped with a liquid-helium cooling pot as 
well as an ion-beam etching system. Ion-beam etching and cryogenic measurements could 
be repeated without breaking the vacuum. The remnant 300-A-thick Au film and Bi-2212 
crystal underneath were further etched, while forming four small stacks on the base mesa. 
The junction size was defined by the width (~ 15 fim) of the base mesa and (~ 12 fim) of 
the electrical pads overlaid perpendicularly to the length of the base mesa. In this etching 
stage we used sufficiently low ion-beam voltage (Heam = 100 V) and current density {heam 
= 0.16 mA/cm^) to enhance the controllability of the etching rate. The height of a small 
stack was controlled by the total etching time tg. In every 1-2 min etching, the chamber 
was pumped out and a stack was characterized in a three-probe-measurement configuration 
[see the inset in Fig. 2(b)]. The Rc(T) was measured using a dc method with a constant 
bias current of 10 /iA. To reduce the effect of external noise, low-pass filters were connected 
to the voltage and the current leads. A recent reportS indicates that, using the similar 
cleaving technique on Bi-2212 crystals, step structures with a height of one half unit cell 
usually forms in the samples in an area of 30 x 30 fim^. We did not observe any evidence 
for the existence of a step in our tunneling results. In this study, we took measurements for 
three different stacks, all showing similar features. We present a typical set of data obtained 
from one of the specimens. 



III. RESULTS AND DISCUSSION 
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A. Rc{T) Curves 



Progressive evolution of Rc{T) by increasing the in-situ etching time up to 30 min is 
presented in Fig. 2(a). Increment of the number of junctions in the stack with increasing 
te was monitored by counting the number of quasiparticle branches in the corresponding 
IV curves. Fig. 2(a) reveals that when we lower the temperature from room temperature 
down to the bulk superconducting transition temperature, the Rc{T) curves show metallic 
behavior at tg < 5 min but gains a slightly semiconducting behavior at te > 6.5 min where 
the value of Rc reaches a maximum at Tc^onset (— 92 K). The value of Rc then starts dropping 
abruptly and reaches a minimum at Tc (~ 86.5 K). Interestingly, for all the values of te, 
Rc remains finite below Tc and gradually increases again up to the secondary resistance 
maximum at 7"^^^^^^^. In the range of T'^^cmset < T < Tc the temperature dependence of Rc 
shows a more pronounced semiconducting behavior. Below T^ onset ^^e value of Rc drops 
sharply, especially for high values of te', and vanishes to 0.1 f2 at T^ (~ 15 K) and below. 

As te increases, T^^^^^^ decreases from ~ 36 K to ~ 17 K, as denoted by the dotted 
guideline in Fig. 2(a). All other characteristic temperatures, Tc^onset, Tc, and T^, however, 
remain almost insensitive to the etching time. The transition width AT^ (~ 0.3~5.0 K) near 
T^ gets sharper with increasing t^ while ATc (~ 2.3 K) near Tc remains unchanged. The 
behavior of Rc{T) near T^ onset ^^^Y sensitive to t^ but is also highly affected by the 

ambient condition. For instance, the peak value of RciT'^^onset) increases rapidly and T'^ is 
suppressed significantly for RF irradiation of power level as low as —30 dBm (not shown). 
Rc{T) near T^'o^^et sensitive to the RF noise that the data taken without filters are 
greatly affected. 

For the sake of clarity we divide the temperature range into three regions: Region I 
(T < T^^ouset), n {n^onset < T < Tc) , and III (T > Tc), as shown in Fig. 2(a). The 
total number of quasiparticle branches obtained from the IV curves in Region I for the 
corresponding etching times are n = 4 {te—6.5 min), 4 (8 min), 5 (10 min), 5 (12 min), 6 
(15 min), 8 (19 min), 9 (22 min), 10 (26 min), and 12 (30 min). 
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The finite value of Rc found in the stack in Region II is known to be always present 
in a three-terminal measurement .iiflEE Since the IJs in the stack below Tc should be in 
zero-resistance state for a low bias current, the finite Rc at temperatures in Region II in this 
three-terminal configuration is attributed to the surface weak junction (WJ) consisting of 
the topmost Cu-0 bilayer in the normal state, and the adjacent inner Cu-0 bilayer in the 
superconducting state.0 In other words, the surface Cu-0 bilayer underneath the normal- 
metal electrode has a suppressed transition temperature, which is essentially T'^ onset- Below 
T^', the surface junction also becomes Josephson coupled and the junction resistance subtracts 
from the value of Rc for a low current biasS In this region, Rc is then essentially the contact 
resistance between the surface Cu-0 bilayer and the Au electrode and its value reduces to as 
small as 0.1 (not shown in detail in the figure), which is at least two orders of magnitude 
smaller than Rc in other temperature ranges. We can thus consider Rc{Tc) as the junction 
resistance of the surface junction. 

We can summarize the above discussion in terms of the tunneling configuration of the 
stack and the corresponding c-axis resistance in the following way. In Region I, the stack 

has Au/D'IDIDID configuration, where D (D') denotes the (suppressed) superconducting 

layer with dx2_y2 symmetry. Rc in this case is essentially the contact resistance Rct between 
Au electrode and D' layer; Rc{T) = Rct ~ 0.1 Q. In Region II, the configuration becomes 

Au/N'IDIDID , where N' denotes the surface layer in its normal state. The corresponding 

resistance is0 

RciT) = Rct + Rn'id{T) ~ Rn'id{T) (1) 



Rn'id{T) 





^ f , d 1 


lim — 


= < lim / 


v-^o dV 


\v-o dV eR'^ J-oo 



-1 



(2) 



where Rn'id{T) is the quasiparticle tunneling resistance of the N'lD surface junction or 
WJ, R'n its junction resistance, and f{E) the Fermi distribution function. The normal- 
ized density of states N{E, A) of the superconducting Cu-0 bilayer with dx2_y2 sym- 
metry can be obtained by averaging over the in-plane angle in k space as N{E, A) = 
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Re[(l/27r) /Q^'^{i?/y^^ — [A(X)s(20)P}c?0]. Here, A is the temperature-dependent energy 
gap of the inner layers. In Region III, the configuration is Au/N'INININ , with the tun- 
neling resistance 



where N denotes the inner Cu-0 bilayer in its normal state, m the total number of IJs (note 
that m = n — 1), and Rn'in (Rnin) stands for the quasiparticle tunneling resistance of each 
N'lN (NIN) junction. 

From Eqs. (P and (||) we note that Rc{Tc) ^ Rn'id{Tc) = R'^ and the value increases 
gradually with increasing tg. The values of Rc{Tc) in Fig. 2(a) are 1.1, 1.2, 1.8, 3.2, 5.2, 
6.0, 7.0, 8.3, 9.0, 9.3, and 9.5 Q for te from up to 30 min [see Fig. 7(a)]. We presume the 
variation of Rc{Tc) is due to shrinking of the junction the etching gets longer, the 

reason for which will be discussed below. 

To study the intrinsic properties of the interlayer coupling or the conduction in the Cu-0 
planes, one needs to eliminate the effect of the variation in the junction area by normalizing 
RciT) with respect to the junction resistance. Rescaled curve Rc(T)/Rc(Tc) for each set 
of curves in Fig. 2(a) corresponding to tg equal to 10 min or longer tends to merge into a 
single curve at temperatures in Region II [see Fig. 2(b)]. The temperature dependence of 
the curve is well in accordance with R^'iniT) expressed in Eq. (H). Note that in Eq. (H) 
only R'^ has te dependence, thus Rn'id/ R^n equivalently RdT) / RdTc) as plotted in Fig. 
2(b) should be independent of t^- The best fit denoted as the dotted curve is obtained for 
Ao = A(0) = 32.6 meV with the assumption of the BCS-type temperature dependence of 



the true BCS-type behavior, for this dx2-y2 symmetry a = 1.45 gives a better fit and is 
consistent with the theoretical prediction.^ The fit turns out very satisfactory in almost all 
temperatures in Region II. The value of Aq is in good agreement with the results of other 
studies, obtained from IV curves of the inner IJs with DID-j unction configurations near the 



Rc{T) = RcT + Rn'in{T) + m X Rnin{T) 



Rn'in{T) +m X Rnin{T) 



(3) 




Although a = 1.74 is valid for 
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liquid-helium temperature,a'Llll2l the scanning tunneling spectroscopy,E3E3 or photo-emission 
spectroscopy,!! using optimally or slightly overdoped samples. 

For comparison, we also draw the best fit curve with an isotropic s-wave gapS with 
Ao = 17 meV (the dashed curve) which gives, however, only a marginal fit at best near 
Tc. Any larger s-wave gap Aq would result in even a poorer fit over the range of Region II. 
Insensitivity of the rescaled Rc{T) to the number of IJs in a stack in Region II indicates that 
the variation of Rc{T), at least for te > 10 min, was caused by the variation in the junction 
area with etching. Note that this insensitivity also implies the number of the N'lD junctions 
does not change with progressive etching; only one N'lD junction exists presumably on the 
surface of the stack. The deviation of Rc{T) from the merging curve for tg < 10 min is due 
to incomplete formation of the stack and will be discussed below. One notices in Fig. 2(b) 
that the curves Rc{T)/Rc{Tc) do not show any merging behavior in Region III, even for long 
etching times te > 10 min. This is consistent with the fact that the Rc{T) in Region III 
include the tunneling resistances of both number of inner IJs and the surface WJ, which are 
scaled by i?„ and R'^ [= RdTc)], respectively. Thus, Rc(T) in Region III is not scaled by 
Rc{Tc) alone. 

The inset of Fig. 3 shows more details of Rc{T) curves in Region III for te = 15, 19, 
22, 26, and 30 min corresponding to n = 6, 8, 9, 10, and 12, respectively. Each curve 
shows a pronounced metallic behavior in the high-temperature region, which may imply 
that the conducting planes are in a highly overdoped regime, gradually turning into a 

semiconducting behavior around T = 170 ~ 180 K with lowering temperatures. Increasing 
te tends to increase Rc{T) in proportion as more IJs are included in the stack. In Fig. 3 we 
plot the single junction contribution to the tunneling resistivity Pc,singie as converted from 
the relation! Rc^singieiT) = [Rc{T,i + j) — Re{T,i)]/j using the geometric parameters, the 
junction area 5* = 12 x 15 /im^ and its thickness d = 12 A, for the number of junctions 
i + j = 8, 9, 12, 12 and i = 6, 6, 8, 6, respectively. Thus, the curves correspond to 
averaging c-axis resistivity over 2, 3, 4, and 6 junctions, respectively. Since Rc(Tc) as seen 
in Fig. 2(a) is stabilized for te ^ 15 min, not much relative error is assumed to be involved 
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in this conversion. Note that the surface WJ contribution is subtracted from pc,singie and 
it contains only the inner IJ contribution. One noticeable feature of pc,singie is that the 
linear-T behavior is significantly diminished in the temperature range of T > 170 K, which 
implies that the linear-T behavior in the inset of Fig. 3 arises mainly from the surface 
WJ. The temperature dependence of Pc,singie for different etching stages merge well into a 
single curve. The dotted curve shows the least fit of the data to the empirical relation 
Pc(T) = (a/T) exp(A*/T) + 6T + c as adopted in Ref. ^ anticipating the occurrence of a 
pseudogap above Tc, with the best fit parameters of a = 146 ± 8 f2cmK, A* = 240 ± 10 
K, b = 0.0250 ± 0.008 Qcm/K, and c = 17.3 ± 1.0 Qcm. The parameter values, especially 



the value of b, in comparison with Fig. 4 of Ref. ^ indicate that, although the linear-T 
dependence is much reduced, the inner stacks are still in a overdoped regime. The surface 
layer, at least with Au protection on it, cannot be less doped than the inner layers. Thus, 
the suppressed superconductivity in the surface layer should be attributed to other causes 
rather than its oxygen loss. The plausible cause is the effect of the proximity contact of the 
surface layer with the normal-metal (Au) electrode, as proposed previously.0 

B. Tunneling IV Curves 

Figure 4 shows the evolution of the quasiparticle branches in the tunneling IV curves 
below as the stack forms with increasing etching time tg- For tg shorter than 5 min, the 
IV curves show anomalous splitting of quasiparticle branches. In the case of te = 5 min 
as in the inset of Fig. 4(a), for example, increasing the bias current causes a sudden drop 
of voltage to the lower- voltage branch near Vdrop — 18 mV, instead of a jump to the next 
higher-voltage branch as usually observed in a high- Tc stack. When the direction of the 
bias change is reversed, voltage jumps around Vjump — 7 mV to a neighboring higher- volt age 
branch. This anomalous branch splitting may have been caused by any remnant Au deposit 
which was not completely removed from the surface of the base mesa during the initial 
stage of the stack formation. The remnant deposit appeared to persist for etching up to 5 



11 



min and partially short the surface Cu-0 bilayer. We believe that the positive temperature 
coefficient in the Rc{T) curves above as well as the deviation of the curves from the single 
merging curve in Region II for < 5 min as seen in Fig. 2(b) was caused by this incomplete 
formation of a stack. 

In Fig. 4(a), for = 6.5 min, one sees four branches develop. As the bias current is 
increased voltage jumps to an adjacent higher- voltage branch at the critical current of each 
branch. Development of this more normal behavior of a stack indicates that the regular 
stack structure started forming at = 6.5 min. Even in this case, the critical current Ic 
(or the return current Ir) for each branch is much different from each other. The IV curve 
for 8-min etching shows similar behavior, while the critical current in each branch reduced 
rapidly (not shown here). The significant discrepancy in (or J^) most likely resulted from 
the differences in the junction area for a series of IJs near the bottom of a stack due to 
inhomogeneous etching around the boundary of a stack. The irregularity in the values of Ic 
and Ir appearing in the high-bias range in the initial etching stages is seen to continuously 
reduce, and both of the quantities approach stable values of their own in the low-bias range as 
te increases [see Figs. 4(b)-(f)]. This behavior implies that the branches with the irregular 
Ic and Ir correspond to incompletely developed junctions near the bottom of the stack. 
When the etching time is increased beyond 10 min, at least a few low-bias branches in the 
IV curves start showing almost the same critical currents {Ic ~ 2.3 ± 0.1 mA) and return 
currents {Ir ~ 70 ± 4 fxA). Thus, a few small stacks are believed to develop fully for te > 10 
min. This conclusion is consistent with what we discussed in relation with the Rc{T) curves. 

If we consider the stack as a c-axis one- dimensional array of the underdamped Josephson 
junctions,@0'ii the McCumber parameterll f3c deduced from the ratio0 Ir/Ic ~ 4/(7? y^) ~ 
0.031 is about 1600, which is about two or three times larger than the results of others, i.e., 
(3c = 500 (Ref . |Tl|) or 300 ~ 700 (Refs. [23|J25| ) . This difference may arise from a high critical 
current density jc — 1300 A/cm^ of our sample. One should note that the first quasiparticle 
branch in the IV curves in three-probe measurements comes from the DTD surface WJ. The 
magnified view of the first branch in a very low-bias range in the inset of Fig. 4(b) shows a 
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hysteretic behavior with a much smaller critical current = 46 ± 3 /iA and return current 

j; = 18 ± 1 /iA- 

The negative dynamic resistance or the "back-bending" for a set of longer etching times 
in Figs. 4(d)-(f) is a generic feature in the IV characteristics of stacked tunneling Joseph- 
son junctions.S Either a nonequilibrium quasiparticle distribution in ultrathin conducting 
layeiili or a Joule heating effeclS has been proposed to explain the feature. Fig. 4 clearly 
demonstrates that the back-bending develops gradually as the number of IJs in the same 
stack increases, attributing it to a high-bias effect. In Fig. 5 we replotted the high-bias 
branches showing the back-bending behavior for the different number of junctions in the 
stack. The back-bending feature is visible, even from the early stages of etching (tg = 15 
min corresponding to n = 6), and becomes more pronounced for longer etching. Remarkable 
feature in this figure is that the voltage-turning points lie well on the 3 mW power-dissipation 
line. Although no symptom of temperature change was detected in the thermometry, this 
strongly suggests that Joule heating could be a likely cause of the back-bending effect. 

Careful examination of Fig. 5 reveals that as te increases a certain quasiparticle branch 
for a constant current bias shift to higher voltage range. This is opposite to what one 
would expect from the nonequilibrium quasiparticle distribution or from the Joule heating 
effect. We infer that it arises from the increase in the junction resistance i?^ and i?„ due to 
shrinking of the junction area with increasing te- 

In Figure 6, we replotted the first low-bias branch of IV curve for tg = 6.5 min which 
was already shown in Fig. 4(a) and the first three low-bias branches of 6 sets of IV curves 
corresponding to the etching times from 10 min to 30 min. We also plotted the anomalous 
branches of IV curves for tg = 3 and 5 min. To eliminate the variation of the junction 
resistance with te, the current axis in each IV curve was rescaled by multiplying Rc{Te) 
obtained from Fig. 2(a). One first notices that in this rescaled plot, with increasing te up 
to 6.5 min, the over-all distribution of branches gradually shifts to higher voltages, and the 
first branch of te = 6.5 min almost coincides with the first stable branch for te = 10 min 
or longer. The sequence of the branch development indicates that the first stable branch 
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arises from the surface WJ in D'lD configuration, and the second and third branches from 
the inner IJs, both in DID configuration. 

If we assume that the D'lD (DID) junction has the same junction resistance R'^ {Rn) as 
the N'lD (NIN) junction, we can rewrite the IV relation for the WJ and IJs as follows. 

1 r°° 

Id'id X < = - / N\E- A')N{E + eV; A)[f{E) - f{E + eV)]dE (4) 
1 r°° 

Idid y<Rn = - N{E; A)N{E + eV; A)[f{E) - f{E + eV)]dE (5) 
e J — oo 

in which A' is the temperature-dependent suppressed energy gap of the surface layer. Both 
i?^ and Rn contain the possible variation of the junction area with etching times. Eqs. @ 
and (P) predict that each IV curve of the WJ or the IJ will show its own scaling behavior 
with respect to different junction resistances R'^ and depending only on the energy gap 
A' and A. The remarkable feature in Fig. 6 is that, when scaled by Rc{Tc) all together, 
all the three branches, for significantly different etching times longer than 10 min, fall into 
three merging curves. This confirms a definite proportionality existing between R'^ and Rn, 
regardless of the etching time tg- The point will become clearer in the discussion in relation 
with Fig. 7(a) below. 

We now determine the value of the zero-temperature energy gap in the surface layer, 
Aq, by fitting the first stable branch to Eq. (^), while adopting Aq = 32.6 meV which was 
already obtained from the fit of Rc{T) in Fig. 2(b) and setting R'^ to be Rc{Tc) obtained 
from Fig. 2(a). The best fit valuell is Aq = 15.3 ±1.1 meV assuming the BCS temperature 
dependence of gap energy with a = 1.45. This value of the surface gap is about one-half 
of the intrinsic gap value for any number of fully developed junctions. From the best fit of 
the second or the third branch to Eq. (^, we obtain i?„ = 1.54 Q for = 30 min, which is 
about one sixth of -R^(30 min)= Rc{Tc] 30 min). It yields the characteristic voltage IcRn ~ 
3.5 mV, which is much smaller than 7rAo/2e = 51.2 mV predicted by the Ambegaokar- 
Baratoff relationil in conventional superconductors or the theoretical prediction Aq/c for 
dx2-y2 symmetry.ii The value i?„ = 1.54 corresponds to p„ ~ 23 i7cm, which is close to 
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the value of Pc,singie ~ 20 — 22 Qcm of Us in the high-temperature region (see Fig. 3). Thus 
the enhanced Rc above i?„ around Tc^onset is beheved to be the quasiparticle contribution 
discussed in Ref. ||. 

Since junctions in a stack switch to the higher resistive state before reaching the gap 
edge of a junction, there are uncertainties in determining both the normal-state resistance 
Rn and the critical current Ic directly from the data. As pointed out in Ref. in a fit to 
Eq. the values of Rn and Aq are interrelated. For instance, one gets smaller fit value of 
Rn for the choice of larger Aq. We eliminated this uncertainty by adopting the measured 
values of Rc(Tc) for R'^, and alternatively determined Aq, Aq, and i?„, respectively, from the 
normalized Rc{T) curves in Region II and IV curves in Region I. This process is possible only 
for repeated progressive measurements in a single piece of crystal with the same material 
parameters Aq and Aq through the entire measurements. 

The consistency of the above fits are further ascertained by cross-checking the etching 
time dependence of the junction parameters obtained by different methods. In Fig. 7(a) 
we plot Rn determined from a fit of the second branch to Eq.(^ for different etching times 
of te = 10 — 30 min. We also plotted Rc{Tc) obtained directly from Fig. 2(a) for different 
tg. Although showing a large difference in the absolute values for a given tg, the two data 
sets obtained by different ways, when normalized by the values at te = 30 min, i.e., Rc{Tc] 
30 min) = -R^(30 min)= 9.52 Q and -Rn(30 min) = 1.54 Q, show an excellent agreement 
with one another for all the ranges of tg values. This fact makes it clearer that a definite 
proportionality exists between i?^ and Rn in any stack-forming etching stages. The value 
of Rc{Tc) starts increasing for tg > 5 min and becomes fully stabilized for tg around 19 
min. The minimal etching time required for the full development of the stack can also be 
examined from the tg dependence of the critical currentsi^ of the WJ and Ig of IJs as 
shown in Fig. 7(b). The critical currents were taken from the first (/^) and the second (Ic) 
quasiparticle branches of the IV curves in Fig. 4. Both values of and Ic are stabilized for 
te around 10 — 15 min. 
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IV. CONCLUSIONS 



In this study, we focused on the evolution of the tunnehng characteristics for stacks con- 
taining a few intrinsic junctions, for temperatures above and below the bulk superconducting 
transition temperature of Bi-2212 single crystals. Three-probe measurements in comparison 
with usual four-probe measurements provided an advantage of enabling us to investigate the 
nature of the surface junction with suppressed coupling strength as well as the conducting 
properties of the inner junctions. Compared with previous attempt by others,! we were 
capable of lowering temperatures even below the Josephson coupling temperature of the 
surface junction using liquid-helium cooling. Our study reveals that Rc{T) curves are scaled 
by the junction resistance of the surface WJ below Tc. On the other hand, the IV charc- 
teristics for both IJs and the WJ in stable etching stage for tf, > 6.5 min are scaled by the 
junction resistance of WJ, which indicates that a definite proportionahty exists between the 
values of junction resistances of the WJ and IJs. 

One of the main findings of this study is that the gap in the surface layer Aq is much 
reduced from the bulk value of IJs, while the junction resistance of the surface WJ (i?^) is 
much larger than that of the IJs (-R„). We obtained Rn from a fit to the second IV branch 
as shown in Fig. 6. Mentioned earlier, in a fit to Eq. an almost similar fit quality can 
be achieved by choosing different values for interrelated parameters Rn and Aq. Instead 
of obtaining R'^ from RdTc), just setting R'^ to be Rn, for example, may be a reasonable 
choice. With this choice, however, the fit gives Aq > 3Ao, i.e., the surface layer has a gap 
about three times larger than the inner layers. Recently it has been claimed that, in the 



underdoped regime, the gap Aq has a larger valueaOt3 even for lower transition temperature 
Tj,'. Thus it may seem plausible to expect widening of the gap at the surface layer while 
Tj,' is suppressed. However, the universal relation between the superconducting gap and the 



transition temperature in Ref. |3^ indicates that the magnitude of the gap for 30 K is at 
best Aq ~ I.SAq, which is far smaller than 3Ao obtained from the fit above. In addition, it 
has also been known that the suppression of the transition temperature of Bi-2212 material 
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itself down to 17 — 36 K cannot be achieved with a high degree of oxygen deficiency, because 
a band with Bi-0 antibonding near the Fermi level acts as a source of holes.ill The picture 
that the surface layer is more underdoped than the inner layers is again contradictory to 
the behavior revealed in Fig. 3, where the opposite was concluded earlier. In this sense, 
we discard the picture of an increased gap with suppressed transition temperature at the 
surface layer, and cling to the original results of fit with a suppressed surface gap. In this 
case, however, we must provide an explanation for the concurrent reduction of and Aq in 
the surface layer, which is apparently in contradiction to the known relationil between the 
two variables. The only way to eliminate the inconsistency is to assume that the concurrent 
suppression of Aq and of the surface layer is not caused by the change in the doping level 



in it, but by its proximity contact to the Au normal electrode as supposed by us in Ref. |TT 
Recently, Manabe et al. observed that, in scanning tunneling spectroscopy measurements 
on Bi-2212, the gap becomes smaller as the tip approaches closer to the crystal surface.!! 
We believe this is consistent with our claim of the proximity-induced suppression of the gap 
on the surface layer in our specimen. Although currently no theory is available to allow 
any quantitative confirmation of our claim, this study apparently strengthens our previous 
arguments. 

Recently there have been tunneling spectroscopic studies on the nature of order- 
parameter symmetry of YBCO-123 and Bi-2212 superconducting materials by employ- 
ing high-Tc superconductor /normal metal /conventional s-wave superconductor hybrid 
junctions.ii'B Studies have reported the existence of significant and finite s-wave compo- 
nents for YBCO-123 and Bi-2212 superconducting materials, respectively. Our study indi- 
cates that in tunneling spectroscopy using the above-mentioned hybrid junctions one should 
be concerned about the change in the superconducting state at the crystal surface, not by 
the oxygen deficiency as one might easily assume, but by the contact of (i^2_y2-wave surface 
layers with a normal layer in a junction. On the other hand, STM tunneling spectroscopy 
may be more affected by the oxygen deficiency at the surface. 

The variation of the junction area along with continuous etching, in fact, acts as a dis- 
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advantage in studying the inherent tunnehng properties in the system. Ironically, however, 
this can be advantageous as in this study, because as illustrated both in Fig. 2 and in Fig. 6, 
we were able to be convinced that, once the data are normalized by proper area-dependent 
parameters like junction resistance, the junctions in any etching stage exhibit no discernible 
variation in their properties, regardless of the number of junctions contained. This leads 
us to the conclusion that not much correlation effect is present among the IJs. Material 
parameters such as Aq and Aq do not change in the whole in-situ etching/measurements 
process, which thus gives stronger conviction of the scaling behavior in our Rc{T) and IV 
data as discussed above. 

The in-situ technique is also extremely useful to study the so-called back-bending effect. 
It convinces one that the effect develops by increasing the number of junctions contained in 
a single piece of crystal, which in turn allows us to draw a power-dissipation line. Compar- 
ing the back-bending feature in the IV branches from different crystals in terms of power 
dissipation would be less meaningful. This technique is expected to provide powerful means 
to examine both the interlayer coupling properties in the superconducting state and the 
in-plane conducting nature in the normal state. 
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FIGURES 

FIG. 1. (a) An optical micrograph of a sample prior to in-situ ion-beam etching. Four small 
stacks were fabricated on the top surface of the base mesa. A-F denote the Au electrodes. Data 
are reported for the configuration where A and D (B and F) were used for biasing currents and 
E and F (A and C) for voltage measurements, (b) Enlarged schematic view of the dashed-circle 
region in (a). 

FIG. 2. (a) Progressive evolution of the c-axis resistance Rc{T) curves with increasing etching 
times; te = 0, 3, 5, 6.5, 10, 12, 15, 19, 22, 26, and 30 min from bottom to top. The contact 
resistance was not subtracted. The characteristic temperatures T^, T'^onseti ^c, and Tc^onset and 
the temperature regions I, II, III are defined in the text. The dotted lines show the tg-dependence 
of Tc and onset- i^) ^ replot of the Rc(T) curves, where Rc{T) was rescaled by Rc{Tc = 86.5 
K). The dotted line is a fit to Eq. (I) for a single normal- metal/insulator /d- wave-superconductor 
(N'lD) surface junction with Aq = 32.6 meV, in comparison with a fit (dashed line) with s-wave 
gap Aq = 17 meV. Inset: a schematic view of the configuration for the in-situ three-terminal 
measurements. 

FIG. 3. The single junction contribution to the tunneling resistivity Pc^singie 

(T) as converted 

from the relation Rc,singie(T) = [RdT, i+j) — Rc(T, i)]/j and the geometry parameters for i+j = 8, 
9, 12, 12 and z = 6, 6, 8, 6, respectively. Inset: the Rc{T) curves above Tc for tg = 15, 19, 22, 26, 
and 30 min, corresponding to n = 6, 8, 9, 10, and 12, the total number of quasiparticle branches 
below T^. 

FIG. 4. Current- voltage characteristics for the etching time (a) te = 6.5 min (n = 4), (b) 10 
min (n = 5),(c) 15 min (n = 6), (d) 19 min (n = 8), (e) 22 min (n = 9), and (f) 30 min (n = 12). 
The substrate temperatures were in the range 13.9 — 16.8 K. Inset of (a): anomalous quasiparticle 
branch splitting for the etching of = 3 min and 5 min. Inset of (b): the enlarged view near zero 
bias which shows a hysteretic behavior of the surface WJ. Arrows indicate voltage changes while 
varying the bias current. 
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FIG. 5. A replot of the high-bias branches showing the back-bending behavior in Fig. 4 
for different numbers of junctions in the stack. The dotted hue represents the 3 mW power hue, 
which passes through the voltage-turning points of the curves, pointing toward the Joule-heating 
mechanism for the back-bending behavior. 

FIG. 6. Replot of the anomalous branches of IV curves for = 3 (dot) and 5 (bullet) minutes 
and the first low-bias branch of IV curve for te = 6.5 (open circle) minutes shown in Fig. 4(a), 
together with the first three low-bias branches of 6 sets of IV curves corresponding to the etching 
times of 10 (solid square), 12 (open square), 19 (solid diamond), 22 (open diamond), 26 (solid 
triangle), and 30 (open triangle) minutes. The dotted line in the first branch (the solid line in 
the second branch) is the fit to the behavior of D'lD surface WJ (a DID inner IJ). The specimen 
temperatures were in the range 13.9 — 16.8 K. 

FIG. 7. (a) The etching time dependence of Rc{Tc) obtained from Fig. 2(a) and the nor- 
mal-state resistance i?„ obtained from the numerical fit to Eq. (^) of the IV curves for the second 
branches in Fig. 6 below T^. For comparison the resistance values were normalized by those for 
te = 30 min; Rc{Tc) = 9.52 $7 and i?„ = 1.54 il. The lines are guides to the eyes, (b) The etching 
time dependences of Ic and at T = 13 — 16 K. The critical currents were taken from the first 
(/^) and the second (Ic) quasiparticle branches of the IV curves below as shown in Fig. 4. The 
lines are guides to the eyes. 



25 



(a) bare surface 

of Bi-2212 base mesa Au electrode 




baked photoresist 



(b) baked 

photoresist 



Au 

electrode thin (< 300 A) 
Au film 




Bi-2212 



base mesa 




5 10 15 20 25 30 
Etching time t (min) 



